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Abstract 

(N 

In the aligned two-Higgs-doublet model, the alignment of Yukawa matrices in flavour space 

OO 

guarantees the absence of tree- level flavour-changing neutral currents, while allowing at 

(N 



> 

X 



the same time for new sources of CP violation, implying potentially large effects in many 
low-energy processes. In this work we study the constraints from exclusive radiative 
B — > Vj decays, where V denotes a light vector meson. The current experimental data 
on the CP-averaged branching ratios and the direct CP and isospin asymmetries are 
analyzed. It is found that, while the branching ratios and direct CP asymmetries do not 
constrain the parameter space much further compared to the inclusive B — > X s d 7 decays, 
complementary constraints can be obtained from the isospin asymmetries A(K*j) and 
A(p7). In addition, correlations between the various observables in exclusive B — > Vj 
and inclusive B — > X s ^ 7 decays are investigated in detail, and predictions are made for 
several so far unmeasured observables. 
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1 Introduction 



Within the Standard Model (SM), flavour-changing neutral current (FCNC) interactions, which 
are absent at tree level, are highly suppressed due to the Glashow-Iliopoulos-Maiani mecha- 
nism thereby providing a very sensitive probe of new physics (NP). In this respect, the 
exclusive and inclusive radiative I?-meson decays induced by the quark-level FCNC transitions 
b — > sj and b — > dj are of particular interest; for a recent review, see e.g. (2]. 

On the experimental side, especially the decays corresponding to b — > s'j transitions are 
known with good accuracy, but branching ratios and even CP and isospin asymmetries have 
been measured for several b — >• d'y decays as well [3j[4]. These results stem mainly from the 
two 5-factory experiments Belle and BaBar; LHCb has however recently started to provide 
data as well [5j[6]. On the theoretical side, while the inclusive decays can be essentially cal- 
culated perturbatively with high precision, the exclusive processes are more complicated due 
to the involved interplay of non-perturbative strong interaction effects (2]. The QCD factor- 
ization (QCDF) approach, which will be adopted in this paper, has provided a systematic 
framework for the treatment of exclusive radiative 5-meson decays [7 10 . For alternative ap- 



proaches, see e.g. 11 15 . Due to these experimental and theoretical improvements achieved in 
recent years, the exclusive b — > s(d) r y decay channels are providing important constraints on 



various NP models 16,17 



Due to its simplicity and being the low-energy limit of some more complete theories, the two- 
Higgs-doublet model (2HDM) has been regarded as a popular extension of the SM since its first 
proposal in the 1970s [18] . In its most general version, the model gives rise to unwanted FCNC 
phenomena due to the non-diagonal couplings of neutral scalars to fermions. Different ways to 



suppress them have been proposed, giving rise to a variety of specific implementations 19 



The tree-level FCNCs can be eliminated requiring the alignment in flavour space of the 



Yukawa matrices coupling to a given right-handed fermion 20 . The aligned two-Higgs-doublet 



model (A2HDM) |2T] results in a very specific structure, with all fermion-scalar interactions 
being proportional to the corresponding fermion masses. The only source of flavour-changing 



phenomena is the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix 22 , appearing in 
the W ± and charged scalar interactions. In addition to the fermion masses and the CKM matrix, 
the Yukawa Lagrangian is fully characterized in terms of three complex parameters qf (f = 



2 



u,d,l), which provide new sources for CP violation. The A2HDM leads to a rich and viable 



phenomenology 21,23,24 , with an interesting hierarchy of FCNC effects, suppressing them in 
light-quark systems while allowing potentially relevant signals in heavy-quark transitions. 

In this paper, we study exclusive radiative I?-meson decays within the A2HDM by employing 
the QCDF approach (7j|8). In addition to branching ratios, we consider the CP and isospin 
asymmetries of these decays. Constraints on the relevant charged-scalar couplings to fermions 
are derived from the current data on these observables. Furthermore, we investigate correlations 
between these exclusive observables and the inclusive B — > X s d 7 decays, and predict several 
observables which have not been measured so far. 

Our paper is organized as follows: In Sec. |2j we first recapitulate the theoretical framework 
for B V'-f transitions and present the physical observables in these decays. In Sec. [3j after 
briefly reviewing the A2HDM, we discuss its contributions to exclusive and inclusive b — > s(d) 7 
decays. In Sec. [4], we present and discuss our numerical results, including the correlations 
between exclusive and inclusive observables, before concluding in Sec. |5j The appendix includes 
a discussion of the relevant input parameters. 

2 B — > V7 decays within the QCDF framework 

In this section, we briefly review the theoretical framework for B — > decays within the 
QCDF approach [7j[8]. For more details, the readers are referred to Refs. [7-10 . 



2.1 The weak effective Hamiltonian 



Within the SM, the weak effective Hamiltonian for radiative b — > (D = d, s) transitions can 
be written as [7] 



a " hf'n^ + M"^] +h.c, (i) 



v/2 

where \[ D ^ = V q bV* D are products of CKM matrix elements. For this result the unitarity 
relation + Ac + X[ D ^ = has been used; we have 



= C^t + ^O^ + J^CiOi, (2) 
= Ci {0\ - Of) + C 2 (£>2 - °t) ■ (3) 
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For b — > s transitions, since the term X^H^ is doubly Cabibbo-suppressed, its contribution 
to the decay amplitude is therefore very small. 

We adopt the operator basis introduced by Chetyrkin, Misiak and Miinz, which is charac- 
terized by a fully anticommuting 75 in dimensional regularization [25] , 

Of = D7„(l - 7s )T> p 7 M (l - 7 5 )T a 6 , 0\ = D^(l - l5 )p pf(l - l5 )b , 

03 = 2^(1-75)6 ^q-fq, O s = 2 ^7^(1 - 75)6 ]T ^W 3 ? , 

<? <? 

4 = 2D 7M (l- 75 )T a & ^T^Tg, O e = 2 D 7mi 7 w 7 m (1 " 7 5 )T a 6 £ ffffl"! . 

r = -g^ m 6 5(7^(1 + 75)6^ , O s = -IjLmtD^T'il + j 5 )bG%, , (4) 

where T a (a = 1,...,8) stands for the SU(3)c generators, mj denotes the 6-quark mass in 
the MS scheme, and e (g s ) is the electromagnetic (strong) coupling constant. The correspond- 
ing Wilson coefficients can be calculated using renormalization-group-improved perturbation 
theory 25 27 . For convenience, we collect in Table [l] their numerical values at the scale 



jj, = 4.4 GeV in leading-logarithmic (LL) and next-to-leading logarithmic (NLL) approxima- 
tion, where Cf = C 7 -C 3 /3-AC 4 /9-20C 5 /3-80C 6 /9 and Cf = C 8 +C 3 -C A /Q+20C 5 -10C 6 /3 
are the so-called "effective coefficients" [25] . 

Table 1: Wilson coefficients at the scale = 4.4 GeV in LL and NLL approximation, using two- loop 
running for a s and the input parameters listed in the appendix. 





C x 


c 2 


c 3 


c 4 


c 5 


c 6 


Cf 




LL 


-0.5166 


1.0263 


-0.0052 


-0.0698 


0.0005 


0.0011 


-0.3126 


-0.1488 


NLL 


-0.3062 


1.0083 


-0.0048 


-0.0840 


0.0003 


0.0009 


-0.3062 


-0.1682 



2.2 Factorization formula for the hadronic matrix elements 

Starting from the weak effective Hamiltonian in Eq. (JT]), the hadronic matrix elements for 
exclusive B — > Vy decays can be written in the heavy 6-quark limit as pip] 

(V(p',eMq,rj)\H®\B(p)) = ^7^(0) [(rf •£*) (p'-q) ~ (v* V) (e*-q) 

(•)) 



where the kinematic variables of initial and final states are indicated in the parentheses, and 
the Bjorken-Drell convention for the Levi-Civita tensor, e 0123 = —1, is adopted. The light 
vector-meson state \V) is defined as 

|p~) , \K*~) for S~-meson decays , 

\V) = I -V2\p°) , y/2\u) , \K*°) for 5°-meson decays , (6) 

\<p) ,\K*°) for 5 s -meson decays . 

With the above convenient parametrization, all the dynamical information is encoded in the 
function 7^(0), the calculation of which constitutes a big challenge (^9)- 

In the QCDF formalism, the function 7^(0) can be computed at leading power in a 
Aqcd/^?> expansion in terms of B — > V transition form factors and light-cone distribution 
amplitudes (LCDAs). Explicitly, the following factorization formula holds [7]: 

7f (0) = T x {V) C f + £ / -^, ± M fdu Mu)T^ ± (u,co), (7) 

l\ c TUB ± J UJ Jq 

where fs and Qb,±{&) denote the B- meson decay constant and LCDAs, and fy and <p±(u) the 
corresponding quantities of the transversely polarized vector meson. The first term is expressed 
in terms of the tensor form factor Ti(0), and corresponds to vertex contributions, where the 
spectator quark in the B meson does not participate in the hard process. The second term, on 
the other hand, incorporates the hard scattering of the spectator quark. Accordingly, they are 
referred to as the "form factor" and the "spectator scattering" term, respectively. 

The hard-scattering kernels C± and T^±(u,uj) in Eq. (7) are perturbatively calculable 
within the QCDF framework. Up to next-to-leading order (NLO) in a s , they have the following 
expansions [7]: 

Cf = + + . . . , (8) 

4:71 

T^ ± (u,u) = T™{u,u) + ^T™{u,u,) + ... , (9) 



where the strong coupling a s should be evaluated at the scale /ib ~ m& in Eq. (|8j) and at /i^ ~ 
(m;,AQCD) 1//2 in Eq. (N), corresponding to the typical virtualities in these two terms. Details 



about the calculation and explicit expressions for the coefficients Cf ,% \ Cj 1 '^ and T±'2(u,u), 
T±'±(u,u) can be found in Ref. 7 . 



In addition to the leading-power contributions given by Eqs. (|8p and it is also well- 
known that specific power corrections, such as the weak annihilation in B — > p7 decays, are 
numerically important [7T[9]. Most importantly, the annihilation topologies have been shown 



to provide the main source of isospin asymmetry in B — > K* r y decays 10 . Despite being 
power-suppressed in Aqco/^b; these terms are still computable within the QCDF framework. 
Thus, they are also included in this paper, and denoted by A7j_ |ann and ATl |hsa- Explicit 
expressions for these terms can be found in Refs. (7 28 . 



However, it should be noted that an endpoint divergence is encountered in the matrix 
element of the chromo-magnetic dipole operator 0$, belonging to the term AT^ |h sa - This is a 
well-analyzed problem for exclusive 5-meson decays within the framework of QCDF and/or its 
field-theoretical formulation, the soft-collinear effective theory (SCET) |29|[30]. Unfortunately, 



there is currently no satisfactory solution to this problem within the QCDF/SCET methods 31 - 



33 . Following the treatment adopted in Refs. 10|28 , we regulate this singularity with an ad-hoc 



cutoff 

rl rl-A h /m B 

du (l + pe ix ) / du, (10) 

Jo Jo 

and take ~ 0.5 GeV, together with < p < 1 and < x < 2tt, to give an estimate of the 
theoretical uncertainty related to this power correction. 

2.3 Observables in B — > decays 

In order to define the physical observables in B — > V'-f decays, it is more convenient to express 

'7 



the decay amplitude in terms of a new quantity Ci , which is defined by |7 



Mi) _ Tt } (°) _ rit^eff , /-.-.x 

67 " Ti(0) " 7 + '" ' 1 } 

where i — t,u refers to the two different CKM factors, and the ellipses denote the subleading 
perturbative and power corrections discussed in the previous subsection. 

In terms of the quantity C? \ the decay rate for B — > V'-f decays can be written as [7J 

T(B^ Vl ) = ^<s(l-^ ^rnlTM 2 \^C? + A? , (12) 

where Gp is the Fermi coupling constant, a era = e 2 /An the fine-structure constant, and S = 1/2 
for p° and u mesons, while S = 1 for the other light vector mesons. Within the SM, the 
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decay rate for the CP-conjugate mode, T(B — > Vj), can be obtained from Eq. (12) with the 
replacement A^ — > A^ *. For b — > s transitions, the dominant contribution comes from 
since the term proportional to A$ is doubly Cabibbo suppressed, whereas for b — > d 
transitions \\^\ is of the same order as |A^|, and the interference between them is the main 
source of CP-violating and isospin-breaking effects. 



With the decay rate in Eq. ( 12 ) at hand, the interesting observables in B — > V'j decays can 



be defined as follows [7[j9]: the CP-averaged branching ratio 
with tb denoting the S-meson lifetime, the direct CP asymmetry 

Acp(B Vl) = TV^-nB^v?) 

V ' T(B — » V7) + T(B — > V7) V ; 

and the isospin asymmetries 

Mir) = r(j?° -> *r-° 7 ) - r(£ + -> k*^) 

1 1} T(B° ^ K*° 1 ) + T{B+ ^ K*+ 1 )' 1 ; 

These observables can be used to test the SM and to probe various NP scenarios. Especially 
the two isospin asymmetries A(fC 7 ) and A(pj) provide complementary information to the 
corresponding inclusive decay modes 16,17]. 



Note that we do not discuss the indirect CP violation. The reason is that this observable 
remains proportional to ma /mi, for a b — > D transition, rendering it very small in the decays 
considered here. Observation of a significant non-zero value would therefore imply NP beyond 
the A2HDM. We note, however, that the available measurements are compatible with zero [4]. 

3 The aligned two- Higgs- doublet model 

3.1 Overview of the A2HDM 

The 2HDM extends the SM with a second scalar doublet of hypercharge Y = |. Thus, in 
addition to the three needed Goldstone bosons, it contains five physical scalars: two charged 
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fields H and three neutral ones ip® = {h, H, A}. The neutral mass-eigenstates p\ are related to 
the original scalar-doublet neutral fields Si through an orthogonal transformation tp® = IZijSj. 

The most generic Yukawa Lagrangian with the SM fermionic content gives rise to FCNCs 
because the Yukawa couplings of the two scalar doublets to fermions cannot be simultane- 
ously diagonalized in flavour space. The non-diagonal neutral couplings can be eliminated by 
requiring the alignment in flavour space of the Yukawa matrices j2l|; i.e., the two Yukawa ma- 
trices coupling to a given type of right-handed fermions are assumed to be proportional to each 
other and can, therefore, be diagonalized simultaneously. The three proportionality parameters 

(/ = u, d, I) are arbitrary complex numbers and introduce new sources of CP violation. 

In terms of the fermion mass-eigenstate fields, the Yukawa interactions of the A2HDM 
read [211 



Cy = H + {u [q d VM d V R - ^ MlVV L ] d + q t vM{P R l) -^^V? $ f M f V R f + h.c. 

(17) 

where V denotes the CKM matrix, V Ri l = are the chirality projectors, and the neutral 
scalar couplings are given by 

yfi = + (^2 + i ^ , yt = n iX + (n l2 - % n i3 ) c . (is) 



From Eq. (17) we can see that, in the A2HDM, all fermionic couplings to scalars are propor- 
tional to the corresponding fermion masses, and the neutral-current interactions are diagonal 
in flavour. The only source of flavour- changing interactions is the CKM mixing matrix in the 
charged-current quark sector. All possible freedom allowed by the alignment conditions is de- 
termined by the three family- universal complex parameters which provide new sources of 
CP violation without tree-level FCNCs (2l|Q 

1 The alignment parameters $f are invariant under global SU(2) transformations of the two scalar doublets, 
4> a <p' a — U a b 0& ; i.e., they are independent of the basis choice adopted in the scalar space. Given an arbitrary 
basis, where the two scalar doublets have vacuum expectation values = (0,v a e l9a ) (a — 1,2), is 

a function of V2/V1, 02 ~~ &i and the (complex) proportionality parameter between the two aligned Yukawa 
matrices coupling to the right-handed fermion [2ll. In the so-called 'Higgs basis' [19], where only one scalar 
doublet acquires a vacuum expectation value, is just the proportionality parameter between the two aligned 



Yukawas coupling to the right-handed fermion Jr 21 
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Higher-order corrections induce a misalignment of the Yukawa matrices, generating small 
FCNC effects suppressed by the corresponding loop factors. However, the flavour symmetries of 
the A2HDM tightly constraint the possible FCNC structures, keeping their effects well below the 



present experimental bounds 20,21 ,23,34]. Although all possible sources of CP violation (from 
the Yukawa matrices as well as the scalar potential) are taken into account in the A2HDM [2l] , 
to the order we are working, the radiative 5-meson decays are only sensitive to charged scalar 
exchange; therefore, the relevant CP-violating effects originate in the parameters <j/ and the 
CKM phase. Thus, in this paper we shall focus only on the phenomenology of the charged-scalar 



Yukawa Lagrangian given by the first term in Eq. (17). 



The presence of flavour-blind phases could induce electric dipole moments (EDMs) at a 



measurable level 35 . Direct one-loop contributions to the light-quark EDMs are strongly 
suppressed by the light quark masses and/or CKM factors. However, this suppression is no 
longer present in some two-loop contributions involving scalar exchanges and a heavy-quark 



loop 36,37. For values of |Im(^*^)| < 1, the predicted neutron and mercury EDMs are 
smaller than the present experimental upper bounds, but they could be within the reach of 
future high-precision measurements. A detailed analysis of EDMs within the A2HDM is in 



progress 



38 



3.2 The A2HDM effects in exclusive B — ► V"f decays 

In a 2HDM without tree-level FCNCs, the NP contribution to b — > s(d) 7 transitions comes only 
from the charged-scalar penguin diagrams. In the approximation of vanishing strange quark 
mass, the resulting effective low-energy operator basis remains the same as in the SM, and 
the charged-scalar effect appears only in the short-distance Wilson coefficients at the matching 
scale nw The charged-scalar contribution to the matching condition has been calculated up 



to NLO independently by several groups 39 43 



Specific to the A2HDM and up to NLO, it is found that the nonzero charged-scalar con- 
tribution to the Wilson coefficients at the matching scale [iw resides only in Cfy 8 (fiw)- F° r 
convenience, we can render explicit their dependence on the couplings s u ,di with the following 
compact form: 



C* s (fiw) — C^ u (fiw) + k«| 2 Ci tUU — s*^d C it ud , (19) 
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where <j*<j d = |<j„||?d|e The short-distance coefficients C itUU and C^ u d are dominated by the 
top-quark contributions, and their explicit expressions at the LO and NLO can be extracted 



from Refs. 39 -43 



From Eq. (19), it can be seen that, depending on the relative phase 0, the combined effect 
of the two terms Cj )MU and Ci tU( i can be rather different. In the calculation of these matching 
coefficients, all the heavy particles (including the top quark, vector bosons and the charged 
Higgs boson) have been integrated out simultaneously at the scale fiw] this is a reasonable 
approximation provided that the charged-scalar mass is of the same order of magnitude 
as m\y and m t . The evolution of the Wilson coefficients from the matching scale fiw down to 
the low-energy scale \x\, ~ rrib remains the same as in the SM. Details about the renormalization 



group evolution and the corresponding solution can be found e.g. in Refs. 25 -27 



4 Numerical results and discussion 

With the theoretical framework presented in the previous sections and the input parameters 
collected in the appendix, we are prepared to present and discuss our numerical results in this 
section. 

4.1 SM predictions and experimental data 

Within the SM, our predictions for the CP-averaged branching ratios as well as the CP and 
isospin asymmetries are collected in Table [2} The theoretical uncertainties are obtained by 
varying the input parameters listed in the appendix within their respective ranges and adding 
them in quadrature, while the experimental data, if not stated otherwise, is taken from the 
Heavy Flavor Averaging Group |4|. For completeness, we also present in Table |2]our predictions 
for the inclusive B — > Xgj and B — > Xdj decays, the theoretical framework of which can 



be found in Refs. 47 -50 for the branching ratios, and in Refs. 51 53 for the direct CP 
asymmetries. 

The main theoretical uncertainties are stemming from the transition form factor Ti(0) for 
the branching ratios, the variation of renormalization scales fit,, fihc for the asymmetries in 
B — >■ fry modes, and the first Gegenbauer moment for the asymmetries in B — > if* 7 
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Table 2: SM predictions for exclusive B — > V7 decays within the QCDF framework. The branching 
ratios are given in units of 10 -6 , the direct CP and isospin asymmetries in units of 10~ 2 . For com- 
pleteness, our predictions for inclusive B — > X s j and B — > Xdj decays are presented as well^. The 
theoretical uncertainties are combined by adding them in quadrature. 



Observable 


n 

Exp. data 4 


SM prediction 


Br(£+ 


42.1 ± 1.8 


40 ill 


Bi(B° K*°-f) 


43.3 ±1.5 


40^ 


A CP (B+ ->• X*+ 7 ) 


18 ±29 


-o.i3i8:2 


A CP (B° -> X* 7 ) 


-0.1 ± 1.5 


fi 00 +0.28 
U.OZ _ .27 


A(K*-f) 


5.2 ±2.6 


4 1 +2 - 4 
-2.0 


Bi(B + -> p+7) 


0.98 t il 


1 ^4+0.53 

J-- ^ -0.48 


Br(5 -+ p° 7 ) 


0.86 ig:a 


0.78 lg£ 


A C p(B+ p+ 7 ) 


-11 ±33 


-9.4l|? 


A CP (S° p° 7 ) 


-44 ±51 


-9.ol|:I 


A(p 7 ) 


-46 +17 

^ u -16 


-7-9 



44 + - 18 

-0.16 



0.57 
-8.6 



-0.30 
-0.25 

-3.0 
-3.6 



*0 



Bt(B s ^K 7) 



*o 



A CP (B S ^K 7) 



1.57: 

-8.4 



-0.60 
-0.53 

+2.9 
-3.4 



Br{B s -> 7 ) 

^cp(^, -> 07) 
Br(£ -> X s7 ) 



33 ±3 



341 ± 22 



44 



-1.2 ±2.8 



51 +16 
01 -14 



308 



h23 
-26 



9 9 +0.6 
4.4 _ 2 g 



Br(5 -+ X, 7 ) 
A CP (5 -> X dT ) 



14.1 ±4.9 



45 



46 



14.2 

-48 



+1.8 
-2.9 

+62 
-14 



t The values given here correspond to a photon-energy cut at E 1 — 1.6 GeV [I]. 



modes. We have added a global 15% uncertainty in all exclusive observables to account for 
non-factorizable effects, not yet included in the QCDF framework. We note that, taking into 
account their respective uncertainties, the predictions for all observables are in good or very 
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good agreement with the data, the only tension appearing in the isospin asymmetry in B — > fry, 
which has however a rather large experimental uncertainty. 

The agreement for the inclusive B — > Xg'j (and B — > X<fy) modes implies very stringent con- 
straints on various NP models [2| [l6|[l~7|[50] . For the direct CP asymmetries, on the other hand, 
due to the hadronic component of the photon, there are quite large theoretical uncertainties, 



lowering the predictive power of these observables 53 . 

In the numerical analysis, we impose the experimental constraints in the following way: each 
point in the A2HDM parameter space corresponds to a theoretical range, constructed as the 
prediction for this observable in that point together with the corresponding theory error. If this 
range has overlap with the 2a range of the measurement, we consider the point allowed. In that 
procedure, the relative theory uncertainty is assumed constant over the parameter space. This 
is a reasonable assumption, since the main theoretical uncertainties are due to the hadronic 
input parameters, common to both the SM and the NP contributions. In order to obtain 
constraints on the charged-scalar Yukawa couplings \q u ,d\i we vary the remaining parameters in 
the ranges m H ± G [80,500] GeV and <p G [0,360°]. 

4.2 B ->> X M 7 decays within the A2HDM 

Since the inclusive B — > Xg^'y decay amplitudes are proportional to the effective Wilson 
coefficient C% s (fib), to facilitate the following discussions, we can decompose it in such a way 
that the dependence on the couplings q u and becomes manifest. In the LO approximation 
and at the scale fib = 2.5 GeV, we have numerically 

CfM - ^ mW {i-0.29& (™^Y + .05W» (20) 

{ V. rn H ± J \ m H ± J J 

from which we can see that, even for comparable \q u \ and |<^|, a stringent constraint on the 
combination q*q<i is expected from these decays. 

Following the discussions in Refs. 23,24 , we show in Figs, [l] and [2] the updated constraints 



on the charged-scalar couplings q u and corresponding to the complex and the real case, 
respectively. We make the following observations: 

• The current experimental data on Br(B — > X^) gives the most stringent constraint on 
the couplings q u and for both the complex and real cases. For ~ tt only a small 
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m H ±[GeV] 
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in H ±[GeV\ 
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Figure 1: Constraints on the charged-scalar couplings $ u and q from inclusive B — > X S) d^ decays, 
plotted in the planes \^d\ ~ 4> (left) and — m H ± (right). The charged-scalar mass m H ± and the 
relative phase (ft are varied in the ranges [80,500] GeV and [0,360°], respectively. 



40 

20 


-20 
-40 





jfcBr(B-»X j7 ) ■ 







-1 








Figure 2: Constraints as in Fig. [TJ but for real couplings plotted in the ? u — plane. 

allowed region remains due to the constructive interference between the SM and the NP 
contributions, similar to the type-II 2HDM, while for <fi ~ the interference between 
them is destructive, and there are two allowed regions, corresponding to relatively small 
NP influence (the lower region) and the case where it is about twice the size of the SM 
contribution (the upper region). 
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This implies for real couplings that while simultaneously large values for real q u and ^ 
with the same signs remain allowed, they are excluded for different signs. 

Taking into account the constraint from the branching ratio, the CP asymmetry does not 
constrain the NP parameter space further. 



As is obvious from Eq. (20), the combination \s*Sd\ is strongly correlated with the charged- 



scalar mass, and large values are only allowed for large m#±. 

Since the branching ratio Br(B — > X s 7) is a key observable, it is interesting to investigate 
its correlations with the other observables. Furthermore, as its effect on the A2HDM is not as 
simple as for the type-II 2HDM, but is merely to strongly correlate the different parameters, 
correlations between observables are significantly affected when imposing the corresponding 
constraint. We show both in Fig. [3j constructed as described above. As the uncertainties of the 
other observables are mostly independent of the one from B — > X s r y 1 the cross for the hadronic 
uncertainties is to be applied to each of these points. Apart from the trivial one between 
Br (B — > X S ( fi)-, we observe mild correlations for most observables, the exception being the two 
direct CP asymmetries. As concluded already in |24] , here an improvement in the experimental 
precision might lead to interesting insights, it would however have to be complemented by 
theoretical progress, given the theoretical error for these observables. 

4.3 B ->> if*7 decays within the A2HDM 

For the exclusive B — > Vj decays, the decay amplitudes are also proportional to the coefficient 



C^d^b) through the quantity df 1 defined in Eq. (11). Accordingly, the general observations for 
the inclusive B — > X s ^ 7 decays apply here as well. However, as we have the isospin asymmetry 
as an additional observable, different constraints on the model parameters are expected. In this 
subsection, we shall discuss B — > K*j decays. 

Imposing the current experimental data on B — > K* / y decays as constraints, we show in 
Figs. [4] and [5] the allowed regions for the charged-scalar couplings, corresponding to the complex 
and the real case, respectively. From these plots, the following observations are made: 

• Constraints on the model parameters from the two CP-averaged branching ratios By(B + — > 
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Figure 3: Correlation plots between the observables in inclusive B — > X S ^"1 decays. The green (dark 
grey) and grey (light grey) regions correspond to the real and complex couplings, respectively. The 
solid and dashed lines denote the experimental central value and the corresponding 2a range, while 
the SM predictions with la error bars are shown by the blue (dark) cross. 

K* +/ y) and Br(5° — > K* 0/ y) are quite similar to those from Br(B — > Xgj), since their 
decay amplitudes are all dominated by the coefficient Cf ff (//&). 

• Since there are still large theoretical and experimental uncertainties for the direct CP 
asymmetries in these decays, again almost no constraints can be obtained from these 
observables, which is the reason why they are not shown here. 



The isospin asymmetry A(iT*7) varies like 1/Cy (/■*&) to first order, and consequently has 
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Figure 4: Constraints from B — > if* 7 decays, plotted in the planes — 4> (left) and — 

m H ± (right). The other captions are the same as in Fig. [TJ 

a different dependence on the relative phase 0, as shown in the third plot of Fig. |4| This 
allows to exclude the large same-sign solutions allowed by the branching ratios. 

• Once constraints from the branching ratios and the isospin asymmetry are combined, the 
allowed parameter space is therefore severely reduced, as shown in the last two plots of 
Fig. [4] for complex, and the last plot of Fig. [5] for real couplings. 

In Fig. [6] the correlations between different observables in B — > K*^y decays are shown. 
The strongest non-trivial ones are between the isospin asymmetry and the branching ratios; 
these are however not very effective within the experimentally allowed range. Very large values 
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Figure 5: Constraints as in Fig. [4| but for real couplings plotted in the q u — ^ plane. 

for A(K*j) correspond to the case where a strong cancellation between the SM and the NP 
contributions to Of occurs, making the remaining parts, such as the annihilation and spectator- 
scattering contributions, relatively important. 

Thus, it is concluded that the data from B — » K*j decays, especially the isospin asymmetry, 
constrain the parameter space in a way complementary to B — > X s 7. A further improved 
measurement of this quantity will therefore be important in constraining NP. 



4.4 B — >• fyy decays within the A2HDM 

For the exclusive B — > /ry decays, since the CKM factors Ai and \[ d ^ are comparable in 
magnitude, both of the two decay amplitudes C? and Q contribute effectively. This feature 
makes these decays particularly interesting in constraining the CKM unitarity triangle and 
searching for physics beyond the SM [8 14 -17 . 



Focusing on the A2HDM effect, we show in Figs. [7] and [8] the current constraints on the cou- 
plings q u and <>d from these decays, corresponding to the complex and the real case, respectively. 
Similarly to the discussions for B — > K*j decays, the correlations between the observables in 
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£> — > p7 decays are shown in Fig. [9] Furthermore, we show in Fig. [10] the correlation between 
the isospin asymmetries in B — > K*^ and B — > /ry. From these plots, we make the following 
observations: 

• From the two branching ratios, currently the charged mode Br(I? + — > p + 7) gives a much 
stronger constraint. The direct CP asymmetries are again found not to be able to put 
any constraints on the model parameters. 

• The isospin asymmetry A(/ry) again shows a different dependence on the relative phase <j) 
than the branching ratios, but also than A(i^*7). This is due to the contribution from the 
extra term proportional to Al , which is associated with a different weak phase arg(V u b). 
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Figure 7: Constraints from B — )■ fry decays, plotted in the planes \<£fid\ ~ $ (left) and \^d\ — 
m H ± (right). The other captions are the same as in Fig. [TJ 

• The combined constraint from the branching ratios and the isospin asymmetry reduces 
significantly the parameter space allowed by the individual observables, showing the im- 
portant role played by A(pry). 

• Even with the constraint from Br (B — > X s '-f) imposed, very large values for the direct 
CP and isospin asymmetries remain allowed, which is again due to the relatively large 
term ~ Al^. This part also includes the formally power-suppressed weak annihilation 
contribution, which is in this case enhanced by the large Wilson coefficients C\p- 

• With the constraint from Br (B — > Xg'j) taken into account, even in the A2HDM the very 
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Figure 8: Constraints as in Fig. [7| but for real couplings plotted in the q u — ^ plane. 

large central value of A(py) cannot be reproduced. Its enhancement would imply smaller 
CP asymmetries in both decay modes. 

• There is a strong correlation between the two isospin asymmetries, and relatively large 
deviations from the SM values remain allowed by the data in the A2HDM, although not 
as large as the present central value of A(pj). 

Thus, similarly to B — y fC*7, the isospin asymmetry A(jry) is a very important observable 
in constraining the charged-scalar Yukawa couplings, and a more precise measurement will yield 
even stronger constraints. A confirmation of the present central value to higher precision would 
challenge the SM as well as most 2HDMs, including the A2HDM. 



4.5 Other B ->> V7 decays within the A2HDM 

— *o 

In this subsection, we discuss the b — >■ d decay modes — > and B s — > K 7, as well as 
the b — > s one B s — > 07. The branching ratio of the latter has very recently been measured 
precisely [5]j6], while for B4 — > the uncertainties remain rather large, and thus only very loose 
constraints on the A2HDM parameters can be obtained. Since neither the CP asymmetries for 
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these modes nor observables for B s — >■ K 7 have been measured so far, we predict them within 
the SM and in the A2HDM. The remaining B — > V'-y modes are predicted to be very tiny within 
the QCDF approach, ~ O(10~ 10 ) [54]. An observed enhancement could imply either NP or 
a breakdown of the method, as these are pure annihilation modes. A quantitative discussion 
does not seem appropriate in this CclSCj clS the A2HDM does not imply large enhancements. 

In Fig. 11 we show the constraints from B s — > 07 and Bd — > wy, where the upper four 
and the lower two plots correspond to the complex and the real couplings, respectively. We 
note that Br(B s — > (fry) is beginning to give competitive constraints compared to the branching 
ratios of the previously discussed modes. 
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Figure 11: Constraints from — > uj~y and B s — > (fry. The upper four and the lower two plots 
correspond to the complex and the real couplings, respectively. 
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From the previous discussions, we already know that Br(5 — > -Xyy), A(K*^y) and A(p^y) 
impose strong, complementary constraints on the A2HDM parameters. Thus, we show in 



Figs. [12p4] the correlations between these three observables and the ones in the modes presently 
discussed, including the ones not yet measured. These can be tested in the near future, given 
the expected experimental progress due to LHCb and the next-generation flavour factories. 
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Figure 12: Correlation plots between ~Bt{B — > X s j) and the observables in Bd — > ury, B s — >• K 7 
and B s — > 07 decays, see also Fig. [3| 

Importantly, for all of the unmeasured CP asymmetries values very different from the SM 
prediction remain allowed. This observation is of special interest for the decay B s — > 07, 
because in this case the SM prediction of a tiny asymmetry is almost unaffected by hadronic 
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uncertainties (see also Ref. (l4j). This fact makes it a key observable not only for the A2HDM, 
but also for every model introducing new weak phases in b — > s transitions. 



5 Conclusions 

The A2HDM, characterized by the alignment of Yukawa matrices in flavour space, guarantees 
the absence of tree-level FCNCs, while introducing at the same time new sources of CP violation 
in the charged-scalar couplings to fermions, implying potentially large effects in many low- 
energy processes. In this paper, employing the QCDF approach, we have studied the exclusive 



24 



I 

O 



1.4 
1.2 
1.0 



BQ, 
CP 



0.6 
0.4 
0.2 
0.0 





1 ■ 

1 














-1 


1 
1 

1 

1 



-80 -60 



-40 -20 

A(py)[%] 



20 



40 




ffl 1.0 



-60 -40 -20 

A(py)[%] 



T 
pq 




-80 -60 -40 -20 

A(py)[%] 



20 



40 




-60 -40 -20 
A(P7)[%] 



20 40 





15 




10 




5 


3 




t 





O 




59, 


-5 






O 


-10 








-15 




10 




5 










t 











-5 


u 










-10 



80 



-60 



-40 -20 

A(py)[%] 



20 



Figure 14: Correlation plots between A (/cry) and the observables in 
B s — » 07 decays, see also Fig. [3} 
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radiative B — >■ V7 decays within this specific NP model. 

With the current experimental data on the CP-averaged branching ratios and the direct 
CP and isospin asymmetries as constraints, we have derived bounds on the charged-scalar cou- 
plings from these decays. It is found that, while the CP- averaged branching ratios cannot give 
additional information to the one extracted from inclusive B — > X s ^ 7 decays, complementary 
constraints on the model parameters can be obtained from the isospin asymmetries A (1^*7) 



and A(p7). Thus, we show in Fig. 15 the final combined constraints from the three observables 
Br (B — > X s j), A(K*^y) and A(/ry), where the first two plots correspond to the complex cou- 
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plings, and the third to real ones. The result confirms explicitly the observation that using the 
exclusive observables can exclude a significant additional part of the parameter space. 
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Figure 15: The combined constraints from Br(B — > X s j), A(K*j) and A(pj), see also Fig. |3j 

In addition, correlations between the various observables in exclusive B — > V7 and inclusive 
B — > X s d 7 decays have been investigated, some of which will become relevant with the advent 
of more precise data. 

Experimental progress for these modes is expected from LHCb and the next-generation 
flavour factories. This will either strengthen the constraints shown here, or will show signs of 
non-standard effects. Of special interest in that respect are the two isospin asymmetries, while 
for the direct CP asymmetries to play a role, also theoretical progress in their calculation would 
be necessary. The exception is the asymmetry in B s — > (jyy, which is cleanly predicted to be 
tiny in the SM and provides therefore important information on any NP model with new weak 
phases in b — > s transitions. 
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Appendix: Input parameters 

In this appendix, we collect the relevant input parameters for our calculation. We restrict 
the discussion to the hadronic quantities which dominate the theoretical uncertainties; basic 
parameters that are known precisely, like meson masses and lifetimes, vector boson masses and 
gauge couplings, can be found in [3]. We use two-loop running for a s throughout this paper. 

CKM matrix elements 

In order to extract values for the CKM parameters A, A, p, and fj, observables insensitive 
to the additional contributions of the A2HDM have to be used. These include the moduli 



of CKM matrix elements from super-allowed (3 decays 55 and semileptonic 5-meson decays 
with light leptons |4|, and the CP-violating angle 7 extracted from tree-dominated S-meson 
decays |56| as more obvious choices. In addition, we use two loop-induced quantities: the 
A2HDM contributions to the ratio Am^/ Am s (2] cancel, and also the indirect CP asymmetry 
in B d — > J/ipKs [4] is negligibly affected. A global fit to these constraints yields 

A — 0.80±0.01±0.01 , A = 0.2254±0.0010, p = 0.14±0.01±0.03 , fj = 0.336±£^±0.010 . 

(21) 
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Heavy quark masses 



For the top-quark mass, the most recent results from the Tevatron and LHC read 57 58 

mj ev = (173.18 ±0.56 ±0.75) GeV , m\ nc = (173.3 ± 0.5 ± 1.3) GeV , (22) 

where both numbers are correlated averages of all measurements performed so far. These results 
still assume the measured value to correspond to the pole mass scheme. For this assumption 
we add an additional theoretical uncertainty of 1 GeV. There are ongoing efforts to determine 
the top mass in a definite scheme; the corresponding uncertainties are not yet competitive, but 
indicate our moderate increase of the uncertainty (compared to the larger difference mf ole — 
m^ s ) to be appropriate. We average the two, which yields finally 

m t = (173.2 ±0.4 ±1.75) GeV. (23) 

For the b- and c-quark running masses in the MS scheme, we take the values from [3] 

fn b (m b ) = 4.18 ± 0.03 GeV , m c (m c ) = 1.275 ± 0.025 GeV . (24) 

To get the corresponding pole and running quark masses at different scales, we use the NLO 



MS-on-shell conversion and running formulae collected e.g. in Ref. 59 . 



Nonperturbative meson parameters 

When discussing exclusive B — > V'j decays, the key quantities are the hadronic parameters 
of the involved mesons, such as heavy-to-light transition form factors and the Gegenbauer 
moments of their LCDAs, as these constitute major sources of uncertainties and have to be 
chosen with the corresponding care. Most of these parameters are currently not directly known 
from experiment and have to be determined by nonperturbative methods like QCD sum rules 
and lattice QCD. 

For the Gegenbauer moments, the relevant values are collected in Table [3] Note that the 
values for the uj meson are not actually calculated, but just assumed to be equal to those of 
the p meson, which is the reason why we doubled the corresponding uncertainties. 



For the heavy-to-light tensor form factors, we use 66 
Fb„^v = f' 



^B q ^V + a l^B q ^V + -jT^B^V + a \r B ^ Yl 
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Table 3: Values for the relevant Gegenbauer moments, extracted using QCD sum rules. 



Meson a lj{l (l GeV) a 1)± (l GeV) a 2 ,||(l GeV) a 2 ,±(l GeV) Ref. 



p 






0.15 ±0.07 


0.14 ±0.06 


60 


61 




K* 


0.03 ±0.02 


0.04 ±0.03 


0.11 ±0.09 


0.10 ±0.08 


62 


64 










0.18 ±0.08 


0.14 ±0.07 


65 



to 0.15 ±0.14 0.14 ±0.12 see text 



which expresses the results of the sum rule calculation in terms of the corresponding decay 
constants and Gegenbauer moments (normalized to their central values in that calculation), 
making the inclusion of updated evaluations of these quantities possible. The coefficients, 
including the remaining uncertainties from the sum rule calculation, can be found in the same 
reference. 

We extract the longitudinal decay constants of the vector mesons, f v = fv, from data. For 
the charged mesons, we can use 



,3 



Br(r- -+ V-v T ) = r T ^- \V UD f /»„ 1 - ^ 1 ± 2^ . (25) 

167T \ m l T J \ m l T J 

Note that this formula holds only in the small width approximation. Regarding the "branching 



ratios" to vector mesons, for K* exists a dedicated analysis 67 , while for p we somewhat 
naively use Br(r — » pu) ~ Br(r — > nixu) — Br(r — > 7r7rz/) non _ res ., using the data from [3]. 

For the neutral vector mesons, the decay constants in question are related to the corre- 
sponding radiative decays. We follow in this extraction mainly 13 .^] However, as we do not 
consider the vector meson mixing in the QCDF analysis, we extract the corresponding decay 
constants with and without taking this mixing into account, using the difference as an addi- 
tional theoretical uncertainty. As the analysis is in addition performed in the isospin limit, 
we do the same for the differences between the decay constants of the charged and neutral p 
mesons. 

The transverse decay constants are not accessible experimentally; they are again typically 
calculated in the frameworks of QCD sum rules or lattice QCD. Generally the calculations 



2 Our value for the extracted decay constants differs from Ref. Il3] as the authors of that paper have got a 
numerical error in the evaluation 68 
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for the ratio with respect to the longitudinal decay constant are more stable, therefore we use 
only the ratios from theory calculations. Note that they are scale dependent, and we take into 
account the (LL) running via = f±(po) (« s (/i) /ots(po))^ 23 • 

For the p and K* mesons, we use the QCD sum rule results from [13]. The <fi meson is 



problematic in these calculations, and we prefer to use the lattice result from 69 . However, 



given the relatively large spread with respect to the calculations in 70-72 , we increase the 
theory error to ~ 3%. The u meson is a special there exists neither a QCD sum rule 

nor a lattice QCD calculation for its decay constant. It is commonly (but somewhat crudely) 
estimated to be equal to the one of the p meson. We use this estimate as well, but increase the 
corresponding uncertainty to account for this. 

In Table [4] we collect the extracted decay constants, as well as the resulting values for the 
transition form factors. In addition to the values given there, we use the ratio /s s //s d = 

to determine the decay constant fs ud - 



1.198 ±0.009 ±0.025 



73 



75 



76 



Table 4: Decay constants and form factors for the relevant decays, for details on the extraction, see 
text. The tensor form factors are given for B u ^ — > p, w, but for B s — > <fi. 



Meson 


/||/GeV 


h(2 GeV)//|, 


T\(0) 


Ref. 


B s 


0.228 ±0.001 ±0.006 






73 


75 


P 


0.215 ±0.001 ±0.006 


0.70 ± 0.04 th 


0.276 ±0.001 ±0.039 


3, 


13 


K* 


0.209 ± 0.007 stat 


0.73 ± 0.04 th 


0.302 ± 0.010 ± 0.051 (B u>d ) 


3. 


13 








0.274 ± 0.009 ± 0.044 (B s ) 








0.229 ±0.002 ±0.002 


0.750 ±0.004 ±0.020 


0.335 ± 0.003 ± 0.043 


3. 


69 




0.188 ±0.002 ±0.010 


0.70 ± 0.10th 


0.237 ±0.003 ±0.055 


1= 





t See text. 
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